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Введение. Статья посвящена разработке математической 
модели уровня подземных вод месторождения цементного 
сырья, расположенного в задонско-елецком водоносном 
горизонте, являющимся основным источником 
хозяйственно-питьевого водоснабжения города Липецка. 
Поэтому на стадии разработки месторождения 
необходимо проводить постоянный мониторинг и иметь 
возможность прогнозирования уровня подземных вод. 
Цель работы — идентификация и исследование 
динамической окрестностной модели с переменными 
иерархическими окрестностями уровня подземных вод, 
позволяющей с достаточной точностью прогнозировать 
значение уровня вод в обследуемых скважинах. 
Материалы и методы. Приведено определение 
динамической окрестностной модели с переменными 
отличающейся 
изменяющимися во времени двухуровневыми 


окрестностными связями между узлами первого и второго 


иерархическими окрестностями, 


уровня. В каждый следующий дискретный момент 
времени узлы окрестностной модели меняют свое 
состояние под воздействием текущих управлений и 
состояний узлов, входящих в их окрестности. В качестве 
частного случая рассмотрена модель с линейными 
функциями пересчета состояний. Параметрическая 
идентификация динамической окрестностной модели 
заключается в нахождении параметров системы для 
каждого узла второго уровня и основана на методе 
наименьших квадратов. 

Результаты — исследования. Разработана линейная 
динамическая окрестностная модель с переменными 
иерархическими окрестностями для прогнозирования 
уровня подземных вод месторождения цементного сырья, 
расположенного в задонско-елецком водоносном горизонте. 
Для параметрической идентификации и моделирования 
функционирования рассматриваемой динамической 
окрестностной модели разработано программное 
обеспечение на языке С++, позволяющее для заданной 
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структуры находить параметры функций пересчета 
состояний узлов, а также прогнозировать поведение модели 
в процессе функционирования. Приведена иерархическая 
структура и проведена параметрическая идентификация 
линейной динамической окрестностной модели уровня 
подземных вод. После выполнения параметрической 
идентификации на обучающей 
математическая модель проверена на контрольной выборке. 
Обсуждение и средние 
относительные ошибки идентификации и прогноза 


выборке данных 


заключения. Полученные 
свидетельствуют об адекватности разработанной модели и 
позволяют рекомендовать ее для прогнозирования уровня 
подземных вод месторождения цементного сырья. 


Ключевые 
месторождение 
окрестностная модель с переменными иерархическими 


слова: уровень 
цементного 


подземных вод, 
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Введение. Задонско-елецкий горизонт является основным продуктивным комплексом известняков, 
пригодным для цементной и металлургической промышленности [1]. Известняки светлого и светло-серого 
цвета, с желтоватым оттенком, средней крепости, трещиноватые, мелкозернистые. Часто встречаются крепкие 
известняки, представленные окремнёнными разностями. Известняки нередко пористые и кавернозные. 
Глинистые породы, пригодные для цементного производства, залегают среди четвертичных отложений, 
покрывающих известняки. 

Вскрышные породы представлены плодородным слоем почвы, суглинками, некондиционной глиной, 
песком и слоем щебенки в верхней разрушенной части известняков. 

Гидрогеологические условия месторождения Сокольско-Ситовское месторождение 
приурочено в геоморфологическом отношении к присклоновой части долины реки Воронеж, чем обусловлена 
гидрогеологическая обстановка в районе месторождения [2]. 


простые. 


Питание водоносного горизонта осуществляется за счет инфильтрации атмосферных осадков и за счет 
поглощения паводкового стока по долинам балок и оврагов. Разгрузка подземных вод происходит в реке 
Воронеж. 

Современное положение уровня подземных вод изучено по режимным наблюдениям в 7 скважинах, 
равномерно расположенных на площади месторождения. Уровень подземных вод в меженный период меняется 
от 102,5 м до 109,7 м. 

В настоящее время задонско-елецкий водоносный горизонт является основным источником 
хозяйственно-питьевого водоснабжения города Липецка. В связи с этим разработка известняков может 
производиться только в необводненной их части с оставлением предохранительного целика не менее 2,0 м, что 
рекомендовано санитарной службой Липецкой области. 

Водоприток в разрабатываемый карьер месторождения возможен только за счет атмосферных осадков. 
В этой связи на стадии разработки месторождения необходимо проводить постоянный мониторинг состояния 
подземных вод задонско-елецкого водоносного горизонта и иметь возможность прогнозирования уровня 
подземных вод. 

В работе для моделирования уровня подземных вод рассматриваемого месторождения применяются 
динамические окрестностные модели с переменными иерархическими окрестностями, позволяющие 
моделировать сложные пространственно-распределенные процессы и объекты, изменяющие свое состояние во 
времени [3-8]. 

Окрестностные модели впервые предложены в конце 90-х гг. ХХ века [9]. В [9-12] приведены 
основные определения и алгоритмы теории окрестностного моделирования. Понятия «окрестности», 
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«окрестностные связи» рассматривается и в [13-14]. В них используются агенты, перемещающиеся по 
окрестностям и взаимодействующие друг с другом по определенным правилам. Сегодня теория окрестностного 
моделирования активно развивается. Появились динамические [6-8], недетерминированные модели [11], 
модели с переменными окрестностями [15]. 

Материалы и методы. В данной работе используются динамические окрестностные модели, 
применяющиеся для прогнозирования состояний пространственно распределенных систем. В них каждый узел 
которых является самостоятельным объектом, функционирующим во времени и связанным окрестностными 
связями с другими объектами системы. Это отличает рассматриваемые модели от широко применяющихся 
сегодня нейронных сетей, которые можно использовать для моделирования процесса функционирования 
каждого объекта или узла в отдельности. Методика окрестностного моделирования предназначена для 
одновременного совместного моделирования и прогнозирования поведения всех элементов распределенной 
системы. 

Далее рассмотрим динамические окрестностные модели «вход-состояние» с переменными 
иерархическими окрестностями. Они отличаются изменяющимися во времени двухуровневыми 
окрестностными связями между узлами первого и второго уровня и линейными функциями пересчета 
состояний. Они позволяют, по сравнению с известными одноуровневыми окрестностными моделями, 
выполнять прогнозирование с более высокой точностью. 

Динамическую окрестностную модель «вход-состояние» с переменными [15] иерархическими 


окрестностями можно задать набором №. = (№ , Г, Ч,Х [0],#) ‚ где: 


1) №=(4,0.,0,,О.„) —щ двухуровневая структура окрестностной модели; даный, — 
множество узлов первого уровня; О, и О, — окрестности связей узлов по состояниям и по управлениям 
соответственно; О„ — иерархические окрестностные связи между узлами. Для каждого узла аеА 


п 
определена своя окрестность по состояниям О [а,]<А и управлениям О[а]<А; О, =чО, [а, ] Ь 
а Е 


п 
О, = ЗО, [а,]. 
Некоторым узлам первого уровня а, е А поставлено в соответствие множество узлов второго уровня 
О. „[а,= {а!,...а*}. Все узлы а, Е О[а,], а, е ОДа, | ‚ оказывают влияние на узлы второго уровня а?’ ЕО. [а]. 
В каждый момент времени определен только единственный активный узел Б=1.....с такой, что 
а’ЕО.„[а,,. 


Узел а’ является активным в момент времени , если для него выполняется заданное условие 


активации /*[!] = ‘тае. Для всех узлов второго уровня а’ О.[а’]=О[а, 1; О[а’]=ОДа,]. 


п 


2) Хе =” — блочный вектор состояний окрестностной модели в текущий момент времени, каждый 
блок которого Х[а,| = Хе А” — вектор состояний в узле а, системы, #=1..... п. 

3) ГЕК" — — блочный вектор управлений в текущий момент времени, каждый блок которого 
Иа]=ИШЙе А” — вектор управлений в узле а, системы, 1=1,...и. 

4) С:Хо хИ, >Х — вектор-функция пересчета состояний окрестностной модели, где Хо, 
множество состояний узлов первого уровня, входящих в окрестность О,; 7, — множество управлений узлов 


первого уровня, входящих в окрестность О,. 
Для узлов первого уровня а, е А функция С. будет иметь вид: 


хи+ьй=ви=У хи+ь]=У 6, (1) 


= 


где 6’ Ао ыН® Г а ХР] — функция пересчета состояний для узла второго уровня а’ : 


Для каждого узла второго уровня а’ еО..[а,,й] функция (6 в текущий момент времени в линейном 


случае имеет вид: 
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ХИ+ЬР]= У &Р,ЛХЬЛ+ У ГР, МИЬ+ 8, (2) 

а, ЕО, [а ] а =О,[аь ] 
где аа, е А (К № ый — Узлы первого уровня модели; ХЬ?]е А” — состояние в узле а’ в момент 
времени # ИЬР]еЕ К” — управление в узле а’ в момент времени & 8’, Ле К”””, э[Р,Ке В" , 


8'[Р]е А” — матрицы-параметры модели. 
Ь Ь > 
Для каждого узла второго уровня а’ ®О.,[а,,Ё] функция С) в текущий момент времени равна нулю, то 


есть Х[+ЬР]= [И =0. 
› Р 

5) [0] ЕР" ° — начальное состояние модели. 

6) { — текущий дискретный момент времени функционирования модели. 

Структура окрестностной модели может быть представлена в виде одного двухуровневого графа 
структуры окрестностной модели — ориентированного графа с двумя видами дуг: по состояниям и 
управляющим воздействиям — или двух ориентированных графов — внешней и внутренней структур. 

Параметрическая идентификация динамической окрестностной модели [16-17] заключается в 
нахождении параметров системы для каждого узла второго уровня и основана на методе наименьших 
квадратов: 


Е=У|ХИ+ЬЙ- 6 па. 


Результаты исследования. Рассмотрим линейную динамическую окрестностную модель уровня 
подземных вод месторождения цементного сырья. Граф внешней структуры модели представлен на рис. 1. Узел 
а’ — внешняя среда. В настоящее время для целей мониторинга уровня подземных вод используются 7 
скважин, которые на рис. 1 соответствуют узлам а,-а,. Управляющие воздействия узлов И[Ё,Й Е Е® состоят 


из количества атмосферных осадков в мм и среднесуточной температуры воздуха в °С за последние 5 суток 
перед измерением уровня воды в скважинах в момент времени Е, 1=1,...,8. Состояния узлов Х|ЬЙеЕ ВР равны 


уровню воды в скважине Г в момент времени #, 1=2,,...,8. 





Рис. 1. Граф внешней структуры окрестностной модели 


Узлы первого уровня окрестностной модели а, (1=2,...,8 ) являются иерархическими: О.„[а,]= {а',а’} 


и соответствуют положительной и отрицательной среднесуточной температуре воздуха. Их структура 
представлена на рис. 2. 
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7Й 





Рис. 2. Граф внутренней структуры узлов первого уровня 


Ь КИ 
Система (2) для каждого узла второго уровня а’ окрестностной модели уровня подземных вод будет 
иметь вид: 
.Б БГ; :6 Ы) БГ; .6 БГ +В 
А а + 8, (3) 
реа Е 0. ЬГ: Е $ = 
где в [1ПеЮ; #[1ПеЕК”; & ПЕК; 1=2.....8; Б=12. 

Для параметрической идентификации и моделирования функционирования рассматриваемой 
динамической окрестностной модели разработана программа на языке С++. Начальными данными для 
идентификации являются структура и обучающая выборка. Программа позволяет находить параметры функций 
пересчета состояний узлов, а также прогнозировать поведение математической модели в процессе 


функционирования. 
После проведения параметрической идентификации на обучающей выборке данных полученная модель 


была проверена на контрольной выборке. Нормированные исходные и модельные значения уровня подземных 
вод для узла а, для обучающей и контрольной выборок приведены соответственно на рис. 3 ирис. 4. 
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Рис. 3. Нормированные исходные и модельные значения уровня подземных вод для узла а, 
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Рис. 4. Нормированные исходные и предсказанные значения уровня подземных вод для узла а, 
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Седых И. А. Прогнозирование уровня подземных вод месторождения цементного сырья на основе динамических окрестностных моделей 
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Средняя относительная ошибка идентификации (прогноза) окрестностной модели вычисляется по 
формуле: 





ма ХИ-+ЬЙ-Х, + 
А= 1 уз = [+ ‚Й 1 + 5Й 100%, 
Мп А,Е+ЪЙ 
где Х„[+ЬЙ — состояние узла а, в т-ом кортеже обучающей (контрольной) выборки; Х„[#+Й — 


модельные значения состояния узла а; М — объем обучающей (контрольной) выборки. 
Средняя относительная ошибка прогноза приведена на рис. 5. 
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Рис. 5. Средняя относительная ошибка прогноза для каждого кортежа контрольной выборки 


Средняя относительная ошибка идентификации составила 0,19%, средняя относительная ошибка 
прогноза — 0,23%, что свидетельствует об адекватности разработанной модели и позволяет рекомендовать ее 
для прогнозирования уровня подземных вод месторождения цементного сырья. 

Обсуждение и заключения. Разработана линейная динамическая окрестностная модель с 
переменными иерархическими окрестностями для прогнозирования уровня подземных вод месторождения 
цементного сырья, расположенного в задонско-елецком водоносном горизонте. 

Для параметрической идентификации и моделирования функционирования рассматриваемой 
динамической окрестностной модели разработано программное обеспечение на языке программирования С+-. 

Реализованная модель с достаточной точностью предсказывает уровень подземных вод в обследуемых 
скважинах и может быть эффективно использована для прогноза уровня подземных вод месторождения 
цементного сырья. 
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